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We  have  developed  a  new  technology  to  predict  capacity  fading  of  lithium-ion  cells  using  the  analysis 
of  their  discharge  curves.  A  cell  capacity  is  deduced  from  the  cell  voltage  window  and  the  cell  discharge 
curve,  which  was  obtained  by  superimposing  discharge  curves  of  the  positive  and  negative  electrodes. 
This  new  technology  consists  of  three  steps:  first,  to  analyze  cell  discharge  curves  measured  during 
the  calendar  life  test  in  order  to  evaluate  deterioration  indicators  such  as  the  positive  and  negative 
usable  capacities  and  the  irreversible  capacities  caused  by  the  side  reactions;  secondly,  to  formulate 
these  indicators  as  the  functions  of  the  test  time;  lastly,  to  produce  the  cell  discharge  curves  by  using 
these  functions  and  deduce  corresponding  cell  capacity,  over  storage  days  for  prediction. 

Prediction  of  the  faded  capacities  of  the  cells  after  a  960  day  calendar  life  test  was  attempted  with  the 
data  up  to  360  days  using  the  technology  described  above  and  the  conventional  method  by  extrapolating 
capacity-fading  data  in  accordance  with  square-root  time  rule.  Consequently,  the  proposed  technology 
predicted  the  faded  capacities  more  precisely  than  the  conventional  method,  particularly  under  severe 
conditions. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  emerging  introduction  of  large-sized  lithium-ion  batteries 
into  various  applications  makes  the  technology  that  can  predict 
their  lifetime  much  more  important  than  before.  The  battery  life 
is  usually  evaluated  by  extrapolating  the  capacity-fading  data  of 
calendar  life  tests  or  cycle  life  tests.  The  extrapolation  method  is 
based  on  the  assumption  that  the  capacity  fading  is  proportional  to 
the  square  root  of  time,  which  is  called  t1/2  rule  [1-4]. 

The  t1/2  rule  was  proposed  by  Peled  as  a  mathematical  model  for 
the  formation  and  growth  of  solid  electrolyte  interphase  (SEI)  [1]. 
SEI  is  formed  by  the  side  reactions  between  intercalated  lithium 
ions  and  the  electrolyte  at  the  negative  electrode  surface  and  con¬ 
sumes  available  lithium  ions.  Therefore,  SEI  growth  is  directly 
related  to  the  capacity  fading,  and  the  t1/2  rule  has  been  applied 
to  the  capacity  fading  of  batteries.  However,  capacity-fading  pre¬ 
dictions  made  by  using  the  t1/2  rule  were  not  always  consistent 
with  actual  data,  especially  in  long  term  or  severe  tests. 

The  inconsistency  between  the  predicted  and  actual  capacity 
fading  seemed  to  be  derived  from  the  assumption  that  capacity 
fading  was  simply  caused  by  SEI  formation,  neglecting  the  effects 
of  other  factors.  Ramadass  et  al.  proposed  to  divide  the  capacity 
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fading  of  the  lithium-ion  cells  into  the  rate  capability  loss,  the  active 
material  loss  and  the  lithium-ion  loss  through  disassembling  the 
deteriorated  lithium-ion  cells  [5].  Though  this  method  enabled  us 
to  decide  active  material  losses  directly,  we  could  not  observe  the 
continuous  deterioration  of  active  material  losses. 

Bloom  et  al.  developed  non-destructive  “differential  voltage 
analysis”  using  the  relation  between  capacity  (Q.)  and  differen¬ 
tial  voltage  (dV/dQ.)  [6].  The  “differential  voltage  analysis”  has  the 
potential  to  precisely  separate  the  discharge  curves  of  the  cells  into 
the  discharge  curves  of  the  positive  and  negative  electrodes,  though 
Bloom  et  al.  have  not  made  use  of  it. 

Zhang  and  White  analyzed  the  discharge  curves  of  lithium-ion 
cells  through  single  particle  model  simulation  and  evaluate  the 
usable  positive  and  negative  active  materials  [7].  Recently,  Smith 
et  al.  have  shown  the  “differential  capacity  analysis”,  which  studies 
discharge  curves  of  dQJdV  vs.  V  [8].  Both  groups  analyzed  dis¬ 
charge  curves  of  the  cells  and  proposed  two  kinds  of  deterioration 
indicators:  the  positive  and  negative  active  material  masses  and 
their  electrode  potentials.  But  both  groups  have  not  determined 
the  quantitative  relationship  between  these  indicators.  Thus,  their 
analysis  results  have  not  been  applied  to  the  capacity  fading  pre¬ 
dictions  by  now. 

We  have  developed  a  mathematically  arranged  “differential 
voltage  analysis”  technology  with  which  we  can  separate  visually 
the  discharge  curves  of  lithium-ion  cells  into  those  of  the  positive 
and  negative  electrodes.  In  the  process  of  separation,  two  kinds 


10142 


K.  Honkura  et  al.  /  Journal  of  Power  Sources  196(2011)  10141-10147 


of  deterioration  indicators,  already  proposed  in  previous  works  as 
shown  above  [7,8],  were  quantitatively  evaluated.  Furthermore,  we 
derived  the  irreversible  lithium-ion  loss  consumed  in  SEI  growth 
from  the  results  of  curve  analysis  [9].  We  treated  three  basic  dete¬ 
rioration  indicators,  irreversible  capacity,  and  the  usable  positive 
and  negative  active  material  losses,  assuming  to  be  independent 
each  other.  Since  this  analysis  is  non-destructive,  we  can  trace  the 
changes  of  these  indicators  during  the  calendar  life  or  cycle  life 
tests.  In  this  work,  we  applied  this  analysis  to  the  calendar  life  test 
results  of  lithium-ion  cells  and  predicted  the  discharge  curves  and 
capacities  of  the  cells. 

2.  Experimental 

2.1.  Calendar  life  tests 

We  fabricated  0.7  Ah  class  1 8,650-type  cylindrical  cells  and  used 
them  to  measure  discharge  curves  during  calendar  life  tests.  The 
active  materials  of  the  cells  were  a  manganese-based  layered- 
structure  oxide  for  the  positive  electrode  and  a  hard  carbon  for 
the  negative  electrode.  The  positive  electrode  contained  7.13  g  of 
active  materials,  while  the  negative  contained  3.22  g.  The  positive 
and  negative  electrodes  were  formed  on  the  both  side  of  aluminum 
and  copper  foils  by  ordinary  slurry  coating  process.  Both  of  the  fin¬ 
ished  electrodes  were  wound  together  with  separators  and  sealed 
into  a  cylindrical  can  with  the  electrolyte  of  1  mol  L-1  lithium  hex- 
afluorophosphate  in  a  mixture  of  organic  carbonates. 

The  initial  discharge  curve  and  initial  discharge  capacities  for 
the  cells  were  measured  between  4.1  and  2.7  V  at  24  mA,  which  was 
regarded  as  small  current  enough  to  obtain  sufficient  capacity  of  the 
cell,  after  initial  conditioning.  The  average  initial  capacity  of  them 
was  0.70  Ah  and  it  was  used  as  the  rated  capacity  by  which  C  rate 
was  defined  during  the  whole  span  of  the  calendar  life  test.  Then 
the  cells  were  charged  up  to  3.63  V  which  is  the  midpoint  of  the 
cell  discharge  curve,  and  then  stored  at  25  °C  or  60  °C.  After  30  days 
of  storage,  they  were  discharged,  charged,  and  discharged  at  less 
than  C/25  at  25  °C  to  measure  the  discharge  curves  and  capacities 
(charge-discharge  test).  They  were  then  charged  up  to  3.63  V  at  1  C 
rate  in  CC/CV  mode  and  stored  again  at  each  specific  temperature. 
The  test  interval  was  changed  from  30  days  at  the  initial  stage  to  60 
days  after  450  days.  This  procedure  was  repeated  until  960  days  of 
storage  so  as  to  prove  the  validity  of  our  capacity-fading  prediction 
technology. 

2.2.  Half-cell  test 

We  measured  half-cell  discharge  curves  of  the  positive  and  neg¬ 
ative  electrodes  at  fresh  state  to  use  them  for  reproducing  the 
discharge  curves  of  the  single  cells.  These  half  cells  consisted  of  a 
pristine  positive  or  negative  electrode  of  1 .8  cm2,  which  were  small 
pieces  of  the  electrodes  with  the  same  specifications  for  18,650 
test  cells,  as  the  working  electrode,  a  separator,  an  electrolyte,  and 
a  lithium  foil  as  the  counter  and  reference  electrodes.  They  were 
assembled  in  an  argon  glove  box  at  room  temperature.  Then  we 
charged  and  discharged  them  at  0.1  mA  (about  C/50  rate)  between 
3.0  and  4.3  V  for  the  positive  half  cell  and  between  0.005  and  1.5  V 
for  the  negative  half  cell. 

Battery  testing  system,  TOSCAT-3000,  Toyo  System  Co.,  was 
used  for  charge-discharge  tests  of  18,650  cells  and  half  cells. 

3.  State  analysis 

3.1.  Discharge  curve  fitting 

The  following  sentences  contain  many  alphabetical  symbols 
with  and  without  suffixes.  Please  consult  the  “List  of  symbols”  sum¬ 


marizing  their  precise  definitions  at  the  end  of  this  paper  when 
necessary. 

The  cell  voltage  VceU  (V)  is  the  difference  between  the  positive 
electrode  potential  Vp  (V)  and  the  negative  electrode  potential  Vn 
(V).  Since  Vceu  depends  on  the  discharged  capacity  Q.  (Ah),  the  dis¬ 
charge  curves  of  QrVceU,  Q-Vp,  and  Q-Vn  will  be  the  functions  of 
Vceii(Q),  VP(Q .),  and  Vn(Q),  respectively.  Note  that  Qis  not  full  capac¬ 
ity  of  cell  Wcen  (Ah),  but  is  discharged  quantity  of  electricity  and 
can  be  obtained  by  just  multiplying  discharge  current  by  discharge 
time.  Supposing  discharge  current  is  small  enough,  we  can  neglect 
the  voltage  drop  due  to  the  internal  resistance  of  the  cell,  and  the 
VCe//(QJ  wiH  be  represented  as  follows: 

UQ)  =  W)-w  a).  (i) 

Vp(Q)  and  Vn(Q.)  were  transformed  into  the  specific  discharge 
curves  of  the  positive  and  negative  electrodes  denoted  as  VsPiP(qp) 
and  VsPin(qn ),  where  q  (Ah  g-1 )  was  the  discharged  capacity  per  unit 
mass  of  the  active  material  and  Vsp  (V)  denoted  the  electrode  poten¬ 
tial.  The  relationship  of  qp  and  qn  to  Q. is  represented  in  Eq.  (2)  with 
the  help  of  other  parameters: 

Q  =  mpqp  —  Sp  =  mnqn  —  8n,  (2) 

where  mp  and  mn  (g)  are  the  usable  masses  of  positive  and  negative 
active  materials  defined  as  the  masses  involved  in  the  discharge 
reactions  and  contributing  to  the  discharge  capacity.  8P  and  8n  (Ah) 
can  be  regarded  as  constant  if  the  following  two  conditions  are 
satisfied  during  discharging  to  measure  Vceu(Q):  (1)  mp  and  mn  do 
not  change,  and  (2)  qp  and  qn  are  varied  not  by  irreversible  side 
reactions  that  cause  mpAqp  and  mnAqn  to  differ,  but  by  reversible 
charge-discharge  reactions  in  which  mpAqp  are  equal  to  mnAqn. 

Fig.  1  illustrates  the  relationship  among  discharge  curves  and 
related  parameters  described  above.  WceU  is  the  full  capacity  of 
the  cell  defined  in  accordance  with  its  own  charge-discharge  con¬ 
ditions.  wp  and  wn  represent  the  full  capacities  per  unit  mass  of 
the  positive  and  negative  active  materials,  which  are  defined  in 
accordance  with  the  capacities  observed  for  the  potential  ranges 
mentioned  in  Section  2.2.  They  are  not  variables  but  constants,  dif¬ 
fering  from  qp  and  qn.  Thus,  mpwp  and  mnwn  are  full  capacities  of 
the  positive  and  the  negative  electrodes  in  the  cell.  8P  and  8n  are 
the  marginal  capacities,  which  remain  uncharged  even  after  the 
lithium-ion  cell  are  charged  to  the  upper  limit  of  the  cell  voltage. 
We  can  also  regard  them  as  the  capacities  that  are  still  discharged 
even  if  the  cell  is  in  fully  charged  state,  Q  =  0. 

Applying  Eq.  (2)  to  Eq.  (1)  transforms  VP(Q)  and  Vn(Q)  into 
Vsp.piQp)  and  Vsp,n{qn\ Then,  Eq.  (1)  will  be 

VCell(Q)  —  Vsp,p{qp)  —  Vsp,n(Qn )•  (3) 
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Fig.  2.  Specific  discharge  curves  of  positive  and  negative  electrodes  and  their  dif¬ 
ferential  forms  measured  with  half  cells. 


Supposing  Vsp,p(qp)  and  Vsp,n(qn)  are  specified,  measured  dis¬ 
charge  curves  of  the  cells  VceP(Q)  can  be  reproduced  using  Eqs.  (2) 
and  (3)  by  adjusting  mp,  mn,  8P,  and  8n  as  free  parameters.  Fig.  2 
shows  examples  of  Vsp,p(qp)  and  Vsp>n(qn)  and  their  differential  form 
dVsPiP(qp)/dqp  and  dVsPtn(qn)ldqn  measured  with  a  half  cell. 

The  reason  why  Vsp{q)  are  transformed  into  the  differential  form 
is  as  follows.  Fitting  of  mp,  mn ,  Sp,  and  8n  in  Eqs.  (2)  and  (3)  some¬ 
times  resulted  in  plural  candidates  that  were  equally  accurate. 
Q-V  discharge  curves  of  the  cells  have  flexibility  and  uncertainty 
to  be  fitted  well  with  wrong  parameters  due  to  their  simple  and 
monotonous  shapes.  If  we  differentiate  both  sides  of  Eq.  (3)  with 
respect  to  Qand  then  replace  dQ with  mpdqp  and  mndqn  on  the  basis 
of  Eq.  (2),  the  differential  discharge  curves  dVceP(Q)/dQwill  be 

dVceU((l)  =  dVp(Q)  dVn(Q) 

dQ  dQ  dQ 

=  J_  dVsp'pfap)  _  J_  dVsp,n{qn )  ,4) 

mp  dqp  mn  dqn 


the  positive  active  material  is  mp[dqp\dt).  mn(dqn/dt)  is  equal  to 
mp{dqpldt)  in  reversible  reactions,  while  mn(dqn/dt)  is  larger  than 
mp(dqpldt)  in  irreversible  reactions  due  to  the  trapped  lithium  ions 
in  SEE  Since  we  assumed  mp(dqp/dt)  as  the  reversible  self  discharge 
rate,  [mn{dqn/dt)  -  mp(dqp/dt)]  corresponds  to  the  irreversible  self 
discharge  rate.  Consequently,  the  irreversible  lithium-ion  loss  dur¬ 
ing  storage  from  t\  to  t2,  A Li/oss  (Ah),  will  be 


h 


(5) 


Although  the  integral  is  difficult  to  do  in  actual  cases  because 
of  the  unknown  functions  m(t)  and  q(t),  we  can  calculate  it  if  the 
decrease  rate  of  m  is  low.  When  this  is  the  case,  replacing  the 
function  m(t)  with  the  intermediate  value  [m(ti)  +  m(t2)]/2  barely 
affects  the  result.  Hence  we  evaluated  the  increment  of  irreversible 
lithium  loss  A Lijoss  from  t\  to  t2  by 


mp(ti)  +  mp(t2) 
2 


(6) 


The  usable  masses  of  active  materials  just  before  and  after 
storage,  m(ti )  and  m(t2),  were  evaluated  by  fitting  dVce//(Q)/dQmea- 
sured  at  t\  and  t2.  The  discharged  capacities  per  unit  mass  of  active 
material  at  t\  and  t2,  q(ti )  and  g(t2),  were  evaluated  in  accordance 
with  the  following  five  steps:  (1)  record  the  open  circuit  voltages 
(OCVs)  of  the  cell  at  t\  and  t2;  (2)  reproduce  dVce»(Q)/dQ  measured 
at  H  and  t2  by  fitting;  (3)  transform  the  fitting  result  of  Q-dV/dQ 
discharge  curves,  i.e.  dVceP(Q)/dQ  dVp(Q)/dQ  dVn(Q)/dQ  into  Q-V 
discharge  curves,  i.e.  Vce//(Q),  VP(Q),  and  Vn(Q);  (4)  estimate  Vp  and 
Vn  corresponding  to  OCVs  using  the  transformed  fitting  results;  and 
(5)  derive  qp  and  qn  corresponding  to  Vp  and  Vn  from  the  specific 
curves  for  V=  Vsp(q). 

4.  Results  and  discussion 


When  the  differential  specific  discharge  curves  of  the  posi¬ 
tive  and  the  negative  electrodes,  dVsPtP(qp)ldqp  and  dVsPin(qn)ldqn, 
are  specified,  the  measured  differential  discharge  curves  of  the 
cells,  dVceP(Q)/dQ  can  be  reproduced  by  Eqs.  (2)  and  (4),  like  the 
case  of  Q-V  discharge  curves  shown  as  Eq.  (3).  dVsp>p{qp)ldqp  and 
dVsp>n(qn)ldqn  have  the  shapes  with  amplified  variation,  peaks,  and 
shoulders  as  shown  in  Fig.  2;  the  values  of  differential  cell  voltage 
dVceuldQ  amplify  as  mp  and  mn  decrease  in  accordance  with  Eq.  (4). 
These  features  provide  us  with  the  clues  to  distinguish  the  most 
reliable  fitting  result  among  some  candidates.  Therefore,  we  can  fit 
dVceP(Q)/dQwith  suitable  mp,  mn,  8P,  and  8n  properly. 

3.2.  Deriving  the  irreversible  lithium-ion  loss  (liioss) 

The  capacity  of  irreversible  lithium-ion  loss  caused  by  the  side 
reactions,  Li/oss,  was  estimated  using  the  reproduced  discharge 
curves  of  the  cells.  The  lithium-ion  cells  were  usually  stored  at 
specific  temperatures  and  periodically  cycled  at  25  °C  during  the 
calendar  life  test.  We  supposed  that  the  side  reactions  consuming 
lithium  ions  to  form  SEI  were  caused  not  by  the  charge-discharge 
test  but  by  the  self  discharge  during  storage. 

When  self  discharge  happens,  some  lithium  ions  come  out  of 
the  negative  active  material  while  some  other  lithium  ions  in  the 
electrolyte  go  into  the  positive  active  material.  The  flow  rate  of 
lithium  ions  from  the  negative  electrode  is  mn{dqn\dt )  (Ah  day-1), 
where  t  is  the  days  for  which  the  lithium-ion  cells  was  stored  at 
a  specific  temperature.  Similarly  the  flow  rate  of  lithium  ions  into 


4A.  Half  cell  test  and  calendar  life  test  results 

The  specific  discharge  curves  and  specific  differential  discharge 
curves  such  as  those  shown  in  Fig.  2  were  measured  by  a  half-cell 
test.  Capacity  changes  in  lithium-ion  cells  during  the  calendar  life 
test  are  shown  in  Fig.  3.  The  horizontal  axis  is  the  square  root  of  the 
storage  days  at  the  specific  temperatures.  As  shown  in  Fig.  3,  the 


Fig.  3.  Capacity  changes  of  lithium-ion  cells  in  calendar  life  test  at  25  °C  or  60  °C  at 
3.63  V. 
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Fig.  4.  Fitting  results  of  differential  discharge  curves  of  lithium-ion  cells:  (a)  in  initial 
state,  and  (b)  stored  at  60  °C  for  720  days. 


change  in  the  capacity  at  25  °C  followed  the  t1/2  rule  represented 
by  the  dotted  line,  while  that  at  60  °C  gradually  deviated  from  the 

t1/2  ruie. 

We  applied  the  state  analysis  we  developed  to  elucidate  the 
reason  for  this  difference.  Fig.  4  shows  the  fitting  results  of  the 
differential  discharge  curves  of  the  cell,  dVce//(Q.)/dQ.  along  with 
dyp(Q.)/d()  and  dVn(Q.)/dQ,  measured  at  25  °C:  Fig.  4a  before  storage 
and  Fig.  4b  after  storage  at  60  °C  for  720  days.  Both  differential  dis¬ 
charge  curves  of  the  cell  were  fitted  to  the  measured  data  very  well. 
Fig.  4a  contains  the  effect  of  irreversible  capacities  caused  by  side 
reactions  in  the  very  first  charge,  which  is  reflected  to  the  large  hor¬ 
izontal  gap  between  the  end  of  the  discharge  curve  of  the  positive 
and  negative  electrode.  When  a  side  reaction  occurs  at  the  negative 
electrode  and  consumes  lithium  ions  coming  from  the  positive  elec¬ 
trode  during  the  first  charge,  the  potential  of  the  negative  electrode 
remains  unchanged,  while  the  potential  of  the  positive  electrode 
rises.  Hence,  the  potential  difference  between  the  positive  and  neg¬ 
ative  electrode  changes,  which  slides  dVp(Q.)/d()  and  dVn(Q.)/dQ.  in 
opposite  direction  horizontally  to  make  the  gap  between  them 
larger.  Comparing  dVce//(Q.)/dQ.  after  720  days  in  Fig.  4(b)  with  that 
before  storage  in  Fig.  4(a),  the  whole  dVce//(Q.)/dQ.  shifted  upward 
while  the  horizontal  width  shrunk,  i.e.  the  differential  voltage  dV/dQ. 
values  increased  as  the  cell  capacity  decreased  from  0.69  Ah  to 
0.52  Ah.  The  upward  shift  is  caused  mainly  by  decrease  of  the  usable 
masses  of  the  positive  and  negative  active  materials,  mp  and  mn. 
The  evaluated  mp  and  mn  were  6.45  g  and  2.02  g  in  the  initial  state, 


Fig.  5.  Utilization  ratios  of  active  materials  in  positive  and  negative  electrodes. 

which  decreased  to  5.77  g  and  1.86  g  after  720  day  storage.  As  mp 
and  mn  decreased,  differential  potential  of  the  positive  and  nega¬ 
tive  electrodes,  dVp/dQ  and  dVn/dQ,  increased  as  shown  in  Eq.  (4). 
Hence,  positive  and  negative  discharge  curves  also  shifted  upward. 
The  horizontal  position  of  dVp(Q.)/dQ.  against  dVn(Q.)/dQ.  changed, 
e.g.  8P  >8n  at  the  initial,  while  8n  >8P  after  720  days.  This  change 
is  presumed  to  be  caused  by  the  side  reactions  at  the  surface  of 
the  negative  active  material  [6].  When  a  side  reaction  occurs  at  the 
negative  electrode  during  storage,  the  intercalated  lithium  ions  are 
consumed  and  the  potential  of  the  negative  electrode  rises,  while 
the  potential  of  the  positive  electrode  remains  unchanged.  It  also 
results  in  increase  of  the  horizontal  gap  between  d\/p(Q.)/dQ.  and 
dVn(Q.)/dQ.  All  the  measured  dVce//(Q.)/dQ.  in  the  calendar  life  tests 
were  well  fitted  in  the  same  way  as  shown  in  Fig.  4.  During  the 
fitting,  every  dVceU(Q)  I  dQ  was  consequently  resolved  into  the  four 
parameters,  mp,  mn,  8P,  and  8n,  which  enabled  us  to  reproduce  the 
original  dVce//(Q)/dQ. 

Fig.  5  plots  the  calculated  utilization  ratios  of  active  materials 
m(t)/m0,  which  are  the  ratios  of  the  mass  of  usable  active  materi¬ 
als  to  the  mounted.  The  utilization  ratios  of  both  the  positive  and 
negative  electrodes  decreased  exponentially  over  the  test  time.  The 
respective  dotted  lines  are  the  fitting  curves  formulated  as  follows 

^=a-0(  l-e-^)  (7) 

m0 

where  m(t)  is  usable  mass  of  the  active  material  as  a  function  of  the 
storage  time  t,  m0  is  the  mounted  mass  of  the  active  material  on  the 
positive  or  negative  electrode,  a  is  the  initial  utilization  ratio,  is 
the  coefficient  of  extent  of  the  decreasing  utilization  ratio,  and  r  is 
the  time  constant  of  the  fading  process.  The  calculated  values  of  the 
fitting  constants  in  Eq.  (7)  are  shown  in  Table  1 .  /3  at  60  °C  was  larger 
than  that  at  25  °C,  and  r  at  60  °C  was  smaller  than  that  at  25  °C. 
This  implies  that  the  usable  active  materials  fade  rapidly  at  high 
temperatures.  The  difference  of  positive  electrode  j3  between  at 
60  °C  and  at  25  °C  is  small,  while  r  at  60  °C  is  around  one  tenth  of  that 
at  25  °C.  This  is  consistent  with  the  quick  decrease  of  m(t)/m0  with 
time  at  higher  temperatures,  for  smaller  r  corresponds  to  quick 
decrease  of  m(t)/mo.  P  for  negative  electrode  at  25  °C  is  a  negative 
value,  which  means  that  m(t)/m0  did  not  decrease  but  increased  as 


Table  1 

Parameters  for  usable  masses  decreasing  in  positive  and  negative  electrodes. 


Electrode 

mol  g 

r/°c 

a 

P 

r/day 

Positive 

7.13 

25 

0.8941 

0.0773 

633.8 

60 

0.9147 

0.1036 

61.1 

Negative 

3.22 

25 

0.6170 

-0.0013 

9.08  xlO5 

60 

0.6155 

0.0745 

1163.5 
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Fig.  6.  Marginal  capacities  of  the  positive  and  negative  electrode. 

time  passed.  Therefore,  the  curve  fitting  for  the  negative  electrode 
at  25  °C  contained  some  errors  caused  by  the  unchanged  values  and 
measurement  errors.  As  for  the  negative  electrode,  r  at  25  °C  is  one 
thousand  times  as  large  as  that  at  60  °C.  This  result  is  consistent 
with  almost  unchanged  m(t)/m0  for  the  negative  electrode  at  25  °C. 

Fig.  6  plots  the  calculated  marginal  capacities  of  the  positive 
and  negative  electrode,  8P  and  8n.  All  dotted  lines  were  the  results 
of  regression  by  power  functions.  As  shown  in  Fig.  6,  8n  increased 
steeply,  while  8P,  decreased  moderately.  To  evaluate  A Li/oss  in 
accordance  with  Eq.  (6),  qp(t )  and  qn(t)  must  be  evaluated  as 
described  in  the  last  part  of  Section  3.2.  Fig.  7a  shows  qp{t)  change 
stored  at  3.63  V  under  25  °C  or  60  °C,  and  Fig.  7b  shows  that  for 


Fig.  7.  Change  of  residual  and  recovered  discharged  capacities:  (a)  positive  elec¬ 
trodes,  and  (b)  negative  electrodes. 


Fig.  8.  Accumulated  irreversible  capacity  losses  by  side  reactions. 

qn(t).  qp  and  qn  of  “residual”  corresponded  to  the  cell  voltage,  which 
were  measured  just  after  storage  of  prescribed  days,  lowered  to  less 
than  3.63  V  by  self  discharge  during  storage.  The  horizontal  axis 
for  “residual”  qp  or  qn  shows  the  days  when  the  previous  storages 
ended  and  the  cell  voltages  were  measured.  On  the  other  hand, 
qp  and  qn  of  “recovered”  meant  that  they  corresponded  to  3.63  V  of 
the  cell  voltage  set  for  the  following  storage.  The  horizontal  axis  for 
“recovered”  qp  or  qn  shows  the  days  on  which  the  cells  were  charged 
to  3.63  V  and  the  following  storages  started.  In  accordance  with  the 
definition,  “residual”  and  “recovered”  discharged  capacities  do  not 
mean  the  quantity  of  the  cell  capacities  to  be  discharged,  but  they 
are  a  kind  of  qp  or  qn  defined  in  Fig.  2.  That  is  to  say,  they  show 
the  distance  from  the  origin  of  horizontal  axis  in  Fig.  2.  Therefore, 
“residual”  included  both  of  reversible  and  irreversible  self  discharge 
capacity,  while  “recovered”  did  not  include  the  reversible.  This  is 
the  reason  why  “residual”  was  larger  than  “recovered”. 

Fig.  8  shows  Liioss(t)  calculated  in  accordance  with  Eq.  (6)  using 
mp(t),  mn(t),  qP(t),  and  qn(t )  shown  in  Figs.  5  and  7,  which  is  the 
accumulated  irreversible  capacity  loss  until  certain  storage  time  t. 
The  dotted  lines  are  fitting  curves  based  on  the  t1/2  rule.  Good  fitting 
of  the  curves  in  the  almost  all  area  in  Fig.  8  indicates  that  Li/oss(t) 
corresponds  to  the  mass  of  lithium  in  the  SEI  film  formed  during 
the  calendar  life  tests,  because  the  t1/2  rule  was  proposed  by  Peled 
assuming  formation  and  growth  of  SEI  [1].  The  non-zero  intersect 
appearing  in  Fig.  8  at  60  °C  may  imply  that  the  SEI  growth  during 
the  first  storage  term  at  60  °C  was  exceptionally  fast. 

4.2.  Capacity-fading  prediction 

All  the  results  and  discussion  on  the  analysis  of  dVceU(Q)ld(l 
described  above  suggest  that  three  indicators  for  the  capacity  fad¬ 
ing  of  lithium-ion  cells,  mp,  mn,  and  Li/oss,  were  the  functions  of  the 
storage  time.  Therefore,  dVce//(QJ/dQ  at  any  storage  time  t  can  be 
reproduced  in  reverse  from  mp(t),  mn(t),  and  Li/oss(t).  The  capacities 
at  any  storage  time  t  can  also  be  calculated  with  Vcell(Q.)  trans¬ 
formed  from  the  reproduced  dVce//(Q.)/dQ. 

We  demonstrated  this  idea  by  the  following  steps:  (1 )  measure 
dVceiiima  UP  t°  360  days  and  analyze  them  into  the  functions 
of  rrip(t),  mn(t),  and  Li/oss(t );  (2)  construct  dVce//(Q.)/dQ.  up  to  960 
days  from  the  functions;  and  (3)  calculate  capacities  of  the  cells, 
Wceiitt),  up  to  960  days  and  compare  the  calculated  and  measured 
capacities. 

Examples  of  the  first  step  have  already  been  shown  above:  mp(t) 
and  mn(t)  decreased  exponentially  in  accordance  with  Eq.  (7),  and 
Liioss(t)  increased  in  accordance  with  the  t1/2  rule  as  shown  in  Fig.  8. 
In  the  second  step,  dVce//(Q.)/dQ.  were  constructed  with  Eqs.  (2)  and 
(4),  which  contains  mp,  mn,  8P ,  and  8n.  Eq.  (2)  is  transformed  as 


10146 


K.  Honkura  et  al.  /  Journal  of  Power  Sources  196(2011)  10141-10147 


follows 


—  (5n  —  8P)  —  -^{mnqn  —  mpqp) 


dt 


/  m  dqn 

dqP 

(m"W 

d Li  loss 

(  dmn 

dt  + 

(qnnr 

dmn 


dmn 


dt  q”  dt 


-dr 


dmp 

dt 


(8) 


The  integral  of  Eq.  (8)  from  t\  to  t2,  supposing  fixed  qp  and  qn  as 
(qp)  and  (qn),  will  be 

A(S„  -  Sp )  =  Li,oss(t2)  -  Liloss(tj )  +  <q„)(mn(f2)  -  m„(t! )) 

-(Qp)(mp(t2)-mp(tJ)),  (9) 


where  Li/oss(0)  was  assumed  to  be  zero  as  the  starting  point. 

The  value  of  ( 8n  -8P)  at  a  certain  storage  time  is  the  sum  of 
A(5n  -8P)  until  the  time  and  the  initial  value  of  (<$n  -8P)  deter¬ 
mined  by  fitting  the  initial  dVceU(Q)ldQ  In  constructing  dVce//(Q.)/dQ, 
8P  was  temporarily  regarded  as  zero.  The  initial  dVce//(Q.)/dQ.  was 
constructed  with  mp( 0),  mn( 0),  and  the  initial  value  of  (<5n  -  8P)  in 
accordance  with  Eqs.  (2)  and  (4).  Then,  the  discharge  curve  of  the 
cells  after  30  days  was  constructed  as  follows.  First,  ( qp )  and  ( qn ) 
in  Eq.  (9)  corresponding  to  the  specific  voltage  of  the  cell  from  0 
to  30  days  were  derived  from  the  constructed  initial  dVce//(Q.)/dQ.in 
the  way  described  in  Section  3.2.  Secondly,  (<$n  -  <$P)  after  30  days 
were  calculated  by  adding  A (8n  -8P)  in  Eq.  (9)  to  the  initial  value 
of  (<$n  -  8P ).  Lastly,  the  discharge  curve  of  the  cells  after  30  days  was 
constructed  with  mp(30),  mn( 30),  and  ( 8n  -8P)  after  30  days.  Dis¬ 
charge  curves  of  the  cells  beyond  30  days  were  constructed  in  the 
same  way. 

Fig.  9  compares  the  predicted  and  the  measured  differential 
discharge  data  stored  for  960  days  at  60  °C.  The  predicted  dif¬ 
ferential  discharge  curves  after  960  days  were  constructed  with 
the  calendar-life  test  data  up  to  360  days  in  the  way  explained 
above.  On  the  other  hand,  the  fitted  ones  were  directly  calcu¬ 
lated  with  the  measured  d\/ce//(Q.)/dQ.  after  960  day  storage  in  the 
way  shown  in  Section  3.1.  Although  the  predicted  and  fitted  pos¬ 
itive  discharge  curves  differed  somewhat,  the  predicted  discharge 
curve  of  the  cell  fitted  the  measured  values  very  well.  The  differ¬ 
ence  between  the  predicted  and  the  fitted  came  from  the  rather 
low  accuracy  of  the  usable  mass  of  the  positive  active  materials. 
Since  the  potential  changes  in  positive  electrodes  were  smaller  than 
those  in  negative  electrodes  within  the  cell  voltage  window,  the  fit- 


Q!  Ah 

Fig.  9.  Predicted  discharge  curves  and  measured  discharge  curves  when  lithium-ion 
cells  were  stored  for  960  days  at  3.63  V  and  at  60  °C. 


Fig.  10.  Measured  capacities  (symbols)  and  the  predicted  by  our  analysis  and  by  t1/2 
rule;  capacities  after  360  days  were  calculated  on  basis  of  test  data  before  360  days. 

ting  results  of  the  positive  electrodes  seemed  to  accompany  greater 
error. 

In  the  third  step,  the  cell  capacities  were  calculated  using  the 
predicted  discharge  curves  of  the  cells  and  the  cell  voltage  win¬ 
dow  prescribed  as  4.09  V  to  2.72  V  for  25  °C  and  as  4.08  V  to  2.73  V 
for  60  °C.  Fig.  10  shows  the  predicted  capacity  curves  with  our 
analysis,  that  with  the  conventional  t1/2  rule,  and  the  measured 
capacities.  Solid  lines  are  predicted  capacities  based  on  our  anal¬ 
ysis,  while  dotted  lines  are  the  predicted  capacities  based  on  the 
t1/2  rule.  Empty  symbols  are  the  capacities  measured  and  used  for 
the  analysis,  and  the  black  ones  are  the  capacities  measured  but 
not  used  for  the  analysis.  All  predicted  capacities  are  calculated  on 
the  basis  of  the  calendar  life  test  data  up  to  360  days.  As  for  the 
results  for  25  °C,  both  the  solid  and  dotted  lines  fit  the  symbols 
up  to  960  days,  which  is  consistent  with  the  almost  constant  mn 
shown  in  Fig.  5.  Although  mp  moderately  decreased  in  this  case  as 
shown  in  Fig.  5,  its  effect  on  the  cell  capacities  is  limited  because 
of  the  small  potential  change  in  the  positive  electrode.  Therefore, 
the  decreases  in  the  cell  capacities  in  this  case  mainly  depend  on 
increases  in  Li/oss  expressed  by  the  t1/2  rule.  On  the  other  hand,  for 
60  °C,  the  dotted  line  based  on  the  t1/2  rule  deviates  from  the  sym¬ 
bols  gradually  after  300  days,  while  the  solid  line  fits  the  symbols 
well  on  the  whole.  In  this  case,  the  decreases  in  mp  and  mn  are 
large  enough  to  be  taken  into  account  as  shown  in  Fig.  5.  Hence  the 
cell  capacities  are  determined  by  the  complex  interaction  among 
decreasing  mp,  decreasing  mn ,  and  increasing  Li/oss.  Therefore,  the 
conventional  t1/2  rule  is  too  simple  to  cope  with  this  complicated 
situation. 

The  new  prediction  method  based  on  the  discharge  curve 
analysis  was  more  accurate  than  the  t1/2  rule  method  when  the 
predictions  were  based  on  the  data  from  zero  to  360  days.  Whether 
this  tendency  stays  the  same  even  if  the  range  of  data  used  for  the 
prediction  is  shortened  or  extended  is  a  significant  question.  The 
answer  is  shown  in  Fig.  1 1.  The  horizontal  axis  shows  the  days  dur¬ 
ing  which  the  data  for  the  analysis  were  measured.  The  vertical  axis 
is  the  deviation  ratio  (%)  to  the  measured  capacity  after  960  days, 
which  is  [Wpred{t)-Wmeas]IWmeas,  where  Wpred{t)  is  the  predicted 
capacity  after  960  days  based  on  the  data  up  to  t  days,  and  Wmeas 
is  the  measured  capacity  after  960  days.  As  far  as  the  25  °C  data  are 
concerned,  both  our  analysis  and  the  t1/2  rule  could  precisely  pre¬ 
dict  the  cell  capacity  measured  after  960  days  if  they  were  based  on 
the  measured  data  up  to  200  days  or  more.  However,  as  for  60  °C, 
the  t1/2  rule  could  not  predict  the  cell  capacity  well  after  960  days 
even  if  the  fitting  range  was  extended  to  840  days.  These  results 
explicitly  show  the  superiority  and  validity  of  the  technology  we 
have  developed. 
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Fig.  11.  Accuracy  of  capacity  prediction  by  our  analysis  and  t1/2  rule;  deviation 
of  measured  capacity  after  960  days  from  predicted  capacity  as  function  of  data 
accumulation  term. 

5.  Conclusion 

We  have  developed  a  new  technology  to  predict  the  capacities  of 
lithium-ion  cells  after  storage  at  25  °C  and  60  °C  using  the  analysis 
of  discharge  curves  of  the  cells.  Periodic  measurement  of  the  dis¬ 
charge  curves  of  the  lithium-ion  cells  during  the  calendar-life  test 
and  their  analysis  enabled  us  to  derive  the  positive  and  negative 
discharge  curves.  The  usable  masses  of  the  positive  and  negative 
active  materials  were  analyzed  by  evaluating  their  discharge  curves 
and  were  shown  to  decrease  not  following  the  t1/2  rule  but  expo¬ 
nentially.  The  exponential  changes  of  usable  masses  revealed  to 
be  much  larger  at  60  °C  than  that  at  25  °C.  The  irreversible  capac¬ 
ities  of  lithium  ions  derived  from  the  analysis  of  discharge  curves 
was  determined  to  mostly  follow  the  t1/2  rule  at  both  25  °C  and 
60  °C.  Consequently,  the  discharge  curves  of  the  cells  at  any  storage 
time  up  to  960  days  were  definitely  reproduced  from  the  data  up 
to  360  days,  by  adjusting  three  deterioration  indicators:  the  pos¬ 
itive  usable  mass,  the  negative  usable  mass,  and  the  irreversible 
lithium  loss.  Our  technology  to  predict  the  capacities  up  to  960 
days  from  the  data  up  to  360  days  reproduced  the  measured  capac¬ 
ities  of  cells  both  at  25  °C  and  60  °C  better  than  the  conventional 

t1/2  ruie. 

Although  all  the  results  of  the  analysis  were  hitherto  consis¬ 
tent,  we  have  to  prove  these  results  qualitatively  and  quantitatively 
by  other  analyses,  such  as  disassembling  the  lithium-ion  cells,  in 
order  to  confirm  the  scheme  formed  by  the  analyses  to  be  true. 
We  applied  our  method  to  merely  one-cell  chemistry  in  this  study. 
Therefore,  one  of  our  important  future  works  will  be  to  prove  its 


availability  for  other  cell  chemistries  in  order  to  generalize  the 
concept  on  which  it  is  based. 

List  of  symbols 

dV(Q)/dQ  differential  discharge  curve 

Liioss  capacity  of  the  irreversible  lithium  loss  in  the  side  reac¬ 
tions,  Ah 

m  usable  mass  of  the  active  materials  contributing  discharge 

capacity,  g 

mo  total  mass  of  the  active  materials  mounted  on  the  elec¬ 
trode,  g 

q  discharged  capacity  per  unit  mass  of  active  material  to  a 

certain  point  of  discharge  process,  Ahg-1 
Q  discharged  capacity  for  the  lithium-ion  cells,  Ah 

t  storage  time,  day 

V  potential  of  the  electrode,  or  cell  voltage,  V 

Vsp  potential  of  the  electrode  per  the  unit  mass  of  active  mate¬ 

rial,  V 

V(Q)  discharge  curve 

w  full  capacity  of  the  electrode  per  the  unit  mass  of  active 

material,  Ahg-1 

W  full  capacity  of  the  lithium-ion  cells,  Ah 
a  utilization  ratio  of  the  active  materials  contributing  to 

capacity  at  the  initial  state 

P  coefficient  of  extent  of  the  decreasing  utilization  ratio  of 

the  active  materials 

8  marginal  capacity  of  the  electrode  remained  uncharged 

in  fully  charged  state  of  the  cells,  Ah 
r  time  constant  of  the  fading  process  of  active  material 

mass,  day 

Subscript 

p  positive  electrode 

n  negative  electrode 

cell  lithium-ion  cell 
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